Nonlinear resistance of 2D electrons in crossed electric and magnetic fields. 
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The longitudinal resistivity of two dimensional (2D) electrons placed in strong magnetic field is 
significantly reduced by applied electric field, an effect which is studied in a broad range of magnetic 
fields B and temperatures T in GaAs quantum wells with high electron density. The data are 
found to be in good agreement with theory, considering the strong nonlinearity of the resistivity as 
result of non-uniform spectral diffusion of the 2D electrons. Inelastic processes limit the diffusion. 
Comparison with the theory yields the inelastic scattering time Tin of the two dimensional electrons. 
In the temperature range T — 2 — lOK for overlapping Landau levels, the inelastic scattering 
rate l/ri„ is found to be proportional to T^, indicating a dominant contribution of the electron- 
electron scattering to the inelastic electron relaxation. In a strong magnetic field, the nonlinear 
resistivity demonstrates scaling behavior, indicating a specific regime of electron heating of well- 
separated Landau levels. In this regime the inelastic scattering rate is found to be proportional to 
T^, suggesting the electron-phonon scattering as the dominant mechanism of the inelastic relaxation. 
At low temperatures and separated Landau levels an additional regime of the inelastic electron 



relaxation is observed: Ti, 



INTRODUCTION 



The nonlinear properties of low-dimensional electron 
systems attract a great deal of attention for its fun- 
damental significance as well as for potentially impor- 
tant applications in nanoelectronics. In response to mi- 
crowave radiation and dc bias, strongly nonlinear elec- 
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states |33l 34 SsJIsgI. 37. 1381 has been reported in 
two-dimensional systems of highly mobile electrons in a 
high magnetic field. There has also been great interest in 
the nonlinear response of quantum ballistic constrictions, 
where the effects of quantum interference, spatial disper- 
sion and electron-electron interaction play essential roles 
[40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 5l|. 

Recent experiments, in which a dc electric field applied 
to highly mobile 2D electrons placed in strong magnetic 
fields, have demonstrated a variety of fascinating non- 
linear phenomena i, S, Ei E m [ii. Oscillations of 
the nonlinear magnetoresistance with a magnetic field, 
which appear at a finite dc bias, have been reported 
[1, & Ell • These interesting oscillations, decaying at 
high temperatures [s^l, are attributed to Landau-Zener 
transitions between Landau levels At substantially 
smaller dc biases another important class of nonlineari- 
ties has been identified [ll|, [H^l • 

In this paper we study in detail the effect of the 
small dc electric field E on the longitudinal resistance of 
two-dimensional electrons in GaAs quantum wells placed 
in a strong magnetic field. In such a magnetic field 
the density of states of the 2D electrons is modulated 
due to the Landau quantization of the electron motion. 



The electric field E decreases the resistance significantly 
0, EH, E3, [ill- The effect, existing in a broad range of 
temperatures, can not be explained by an increase of the 
electron temperature due to the heating by the electric 
field E 
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54| . In the paper [ll| the effect is attributed 
to a non-uniform spectral diffusion of the 2D electrons 
induced by the electric field The spectral diffusion 
produces a specific distribution of 2D electrons in the 
quantized spectrum, which is significantly different from 
the canonical Fermi-Dirac form. In fact the observed 
strong nonlinearity is result of the deviations of the elec- 
tron distribution from the Fermi-Dirac function. The 
effect is considerably enhanced in electron systems with 
high mobility and high electron density. The high elec- 
tron mobility provides strong absolute variations of the 
density of states and the spectral diffusion with electron 
energy, increasing appreciably the magnitude of the non- 
temperature deviations. The high electron density pro- 
vides substantial decrease of the electron-electron scat- 
tering, which makes the relaxation of the deviations to 
be weak. 

Effects of an electric field E on the resistance of two di- 
mensional electrons placed in str ong magnetic fields have 
been studied in many works 5^, 56 1. Substantial part of 
these studies was focused on an effect of the electric field 
E on an amplitude of quantum oscillations of the resistiv- 
ity. The quantum (Shubnikov de Haas, SdH) oscillations 
are result of the quantization of the electron spectrum in 
strong magnetic field [H^l ■ The amplitude of the oscilla- 
tions depends significantly on the electron temperature 
571 |58| . It has been found that the amplitude of the 
SdH oscillations decreases with the electric field E [sst - 
The effect is attributed to an increase of the electron 
temperature due to the electric heating. The expla- 
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nation is based on an assumption that the surplus of the 
Joule energy provided by the electric field E is rapidly 
shared among the carriers through electron-electron in- 
teraction, establishing the thermal (Fermi-Dirac) distri- 
bution at an elevated temperature Tg (sol . |6C| . The 
Tg approximation works well in systems with a strong 
electron-electron scattering. It ignores any deviations of 
the non-equilibrium electron distribution from the Fermi- 
Dirac form. The approximation has been widely and suc- 
cessfully used for 2D electron systems with low electron 
density and/or mobility [H^. We note, however, that 
a substantial discrepancy between the temperature Te, 
obtained from the analysis of the amplitude of the quan- 
tum oscillations in the Tg approximation, and the one 
obtained, using another experimental method, has been 
reported in GaAs 2D systems with a high electron mo- 
bility |60i]. 

Despite the apparent applicability of the Tg approxi- 
mation to the overheated electron systems, recent studies 
have revealed an inadequacy of the temperature descrip- 
tion of the nonlinear transport of highly mobile 2D car- 
riers [H Hi, 113 . Instead of the Te approximation in this 
paper we use a different approach [28|. Below we eval- 
uate the distribution function, using an equation of the 
spectral diffusion. In the computations any assumptions 
regarding the shape of the electron distribution function 
are relaxed. In contrast to the approximation the 
new approach to the heating via the direct evaluation of 
the electron distribution function is more universal and 
accurate. It takes into account, in principle, both the 
broadening ("temperature" increase) of the distribution 
function and the deviations of the distribution function 
from Fermi-Dirac form in response to the electric field 
E. The later appears to be the dominant source of the 
strong nonlinearity observed in highly mobile 2D electron 
systems at small electric fields. 

The spectral diffusion is limited by an electron inelas- 
tic relaxation, which moves the electron system back to 
thermal equilibrium. It opens new possibilities to study 
inelastic processes and nonlinear electron kinetics of low 
dimensional systems. In the present paper we explore 
these possibilities. We study the effect of electric fields 
on the resistivity in a broad range of magnetic fields and 
temperatures. We compare the experimental results with 
numerical simulations of the spectral diffusion. The com- 
parison gives the inelastic scattering time of 2D electrons 
in a broad range of magnetic fields and temperatures. 

In the temperature interval T = 2 — IQK for over- 
lapping Landau levels, the inelastic scattering rate 1/Ti„ 
is found to be proportional to the square of the tem- 
perature, indicating the dominant contribution of the 
electron-electron interaction into the relaxation of the 
electron distribution function. At a strong magnetic field, 
at which Landau levels are well separated, the nonlinear 
resistance demonstrates an interesting scaling behavior. 
In this regime at high temperatures the inelastic scat- 



tering rate is found to be proportional to T^, indicating 
leading contribution of the electron-phonon scattering to 
the inelastic relaxation. At low temperature and sepa- 
rated Landau levels an additional regime of the inelastic 
electron relaxation is observed: ^ 

The paper has the following organization. The " Exper- 
imental Setup" section presents the main kinetic parame- 
ters of samples and details of the experiment. The "The- 
ory and Numerical Simulations" section presents basic 
components of the theory and discusses essential steps 
used to calculate the longitudinal resistance. Experi- 
mental results and a comparison with numerical simu- 
lations are presented in the section " Results and Discus- 
sion" . Section " Conclusion" contains a summary of the 
research. 



EXPERIMENTAL SETUP 

Our samples are high-mobility GaAs quantum wells 
grown by molecular beam epitaxy on semi-insulating 
(001) GaAs substrates. The width of the GaAs quan- 
tum well is 13 nm. Two AlAs/GaAs type-II superlattices 
grown on both sides of the well served as barriers, pro- 
viding a high mobility of_2D electrons inside the well at 
a high electron densitv|62^. Two samples (Nl and N2) 
were studied with electron density ni — 12.2 x 10^^ m~^, 
712=8.2 xl0i5 (m-2) and mobility ni= 93 mVVs, ^2=85 
(m^/Vs) at T=2.7K. At higher densities the cyclotron ra- 
dius re of 2D electrons at Fermi level is larger. As it is 
shown below, this increases the spectral diffusion and the 
nonlinear response in strong magnetic fields. 

Measurements were carried out between T=0.3K and 
T=30K in a He-3 insert in a superconducting solenoid. 
Samples and a calibrated thermometer were mounted on 
a cold cooper finger in vacuum. Magnetic fields up to 1 
T were applied perpendicular to the 2D electron layers 
patterned in a form of d=50 fj,m wide Hall bars with 
a distance of 250 fim along the bars between poten- 
tial contacts. A schematic view of experimental setup 
is shown in Fig[TJ To measure the resistance we have 
used the four probes method. Direct electric current Idc 
(dc bias) is applied simultaneously with an ac excita- 
tion lac through the same current contacts (x-direction) . 
The current contacts are placed far away from the mea- 
sured area at a distance of 500 /im, which is much greater 
than the inelastic relaxation length of the 2D electrons 
Lin = {DTiny/"^ ^1 — 5 fim (see below). The later in- 
sures that possible nonlinearities near the current leads 
provide negligibly small contribution to the total nonlin- 
ear response measured in the experiments. 

Experiments are done at fixed magnetic fields corre- 
sponding to maximums of the Shubnikov de Haas oscil- 
lations. At this condition the Fermi level is located at 
a maximum of the density of states and contributions of 
the edge states to the total electron transport is small. 
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FIG. 1: Schematic view of experimental setup. Studied 2D 
electron system is etched in the shape of a Hall bar. White 
area schematically presents the details of the Hall bar: the 
width and the length of the measured part of the sample are 
d =50 /xm and L =250 /im. Direct current Idc is applied 
simultaneously with ac current lac through current contacts 
formed in the 2D electron layer. The longitudinal ac voltage 
Vac is measured between potential contacts displaced 250 fim 
along each side of the sample. 

Below we consider the density of the electrical current 
across the samples to be a constant. 

The longitudinal voltage Vac was measured between 
potential contacts (displaced along the x-direction) us- 
ing a lockin amplifier with 10 Mfi input impedance. In 
the experiments the potential contacts provided insignif- 
icant contribution to the overall nonlinear response due 
to small values of the contact resistance (about Ikil) and 
negligibly small electric current flowing through the con- 
tacts (< 0.1 nA). 

The differential longitudinal resistance r^x = Vac/ lac 
is measured at a frequency of 77 Hz in the linear regime. 
In the experiment a dependence of differential resistance 
Txx — dVxx/dl on the dc bias Idc is measured. The re- 
sistance Rxx of the sample is obtained by an integration 
of the differential resistance: Rxx = {J fxxdl)/ Idc- In 
the paper we compare the resistance Rxx with numerical 
calculations based on recent theory [28| . 

Experiments are done in a classically strong magnetic 



fields {iVcTtr ^ 1), where the uJc is cyclotron frequency 
and Ttr is the transport scattering time. At this condition 
the electric current density J = {Jx, 0) directed along the 
x-axes is almost perpendicular to the total electric field 
E = {Ex,Ey), where Ex Ey The magnitude 

of the Hall electric field Eh = Ey directed along the y- 
axes is almost equal to the magnitude of the total electric 
field \E\. Below we consider the magnitude of the Hall 
electric field Eh to be equal to the magnitude of the 
total electric field E applied to the samples. The local 
Joule heat injected into the 2D systems per second can be 
evaluated with an accuracy better than 2% as: Jx- Ex = 
(o-xxEx + cFxyEy) ■ {axx/crxy)Ey « axx ■ Ejj, where a is 
the conductivity in the strong magnetic field. 

In our experiments the Hall voltage Vxy is recorded 
simultaneously with the longitudinal voltage Vxx- Ob- 
served variations of the Hall conductivity Uxy and the 
Hall electric field Eh with the dc bias were below 1%. 
These variations yield a negligibly small contribution to 
the overall dependence of the longitudinal conductivity 
Uxx on the dc bias. This contribution are ignored in the 
comparison between the experiment and the theory. 

THEORY AND NUMERICAL SIMULATIONS 

In this section we present basic parts of the theory [2^ 
and details of the numerical calculations of the nonlin- 
ear resistivity. The theory considers nonlinear electron 
transport in a strong magnetic field. In the magnetic 
field the electron spectrum is quantized and the density 
of states oscillates with the energy. The period of the os- 
cillations is the cyclotron energy TilOc- The width of the 
Landau levels is F = h/Tq, where Tq is quantum scatter- 
ing time. At low temperatures the time Tq is determined 
by an elastic impurity scattering of the 2D electrons. At 
small quantized magnetic fields the electron spin splitting 
is much smaller the level width F . The spin splitting 
is neglected in the paper. 

The net longitudinal conductivity of the 2D electrons 
fn; — (J XX is a sum of conductivities cr(e) of the levels with 
energy e over all possible energies, weighted with the first 
derivative of the distribution function df jdt [H^l : 

- y" a(e)(-9//ae)de, (1) 

In the leading approximation for a classically strong 
magnetic field the longitudinal conductivity cr(e) at an 
energy e reads [2^ : 

a(e)=az,^2(e), (2) 

where ub = e^v^v^pj^uy^Ttr is the dc Drude conductivity 
in a strong magnetic field _B, v(t) = i'(e)/vQ is dimen- 
sionless density of states (DOS), Ttr and vq = m/irTi' are 
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transport scattering time and the density of states at zero 
magnetic field and vj? is tfie Fermi velocity. The approxi- 
mation neglects effects of the electric field on the electron- 
impurity collision, which yields a negligibly small correc- 
tion to the nonlinear resistance at small electric fields ^28], 
The dominant nonlinear effect is due to a non-trivial en- 
ergy dependence of the distribution function /(e), which 
is a result of non-uniform spectral diffusion of the 2D 
electrons in response to the total dc electric field E ap- 
plied to the system. 

Due to conservation of total electron energy eg in 
the presence of the external electric field E and the 
elastic electron-impurity scattering, the kinetic energy 
of an electron ex depends on the electron position r: 
^k{t) — ~ eEr. As a result of the energy conserva- 
tion, the diffusion motion of the electron in real space 
originates a diffusion of the electron kinetic energy in the 
energy space. The diffusion generates a spectral electron 
flow from occupied electron levels below the Fermi energy 
to empty states above it. The coefficient of the spectral 
diffusion D^{e) is proportional to the coefficient of the 
spatial diffusion D{e) = Vpi>{e)/2uj'^Ttr = r^v{€) /2Ttr'- 
D,{e) = {eEYD{e) - {5rf. The spectral diffusion is 
proportional to square of the cyclotron radius rc and 
the normalized density of states The spectral diffu- 
sion is most effective in the center of the Landau levels, 
where the density of states is high, gradually decreases 
away from the center and is suppressed considerably be- 
tween Landau levels, where the density of states is small. 

The spectral diffusion is described by the Fokker-Plank 
type equation [ist : 



'dc 



^ + E' 
dt t'ot^(e) 



/(e) "/tW 



(3) 



The left side of the equation describes the spectral 
diffusion of a spherical part of the electron distribution 
function / induced by the electric field E in the pres- 
ence of the elastic impurity scattering. The higher angu- 
lar harmonics of the distribution function provide much 
smaller contributions to the net function /, due to much 
faster temporal relaxation. These are neglected in the 
eqini The right side of the equation describes the in- 
elastic relaxation of the distribution function toward the 
thermal equilibrium expressed by Fermi-Dirac function 
fri^)- The inelastic relaxation is taken in, so-called, r 
approximation of the inelastic collision integral. Validity 
of the approximation is supported theoretically in the 
high temperature limit kT ^ Uuc [1^. Below, in the 
numerical calculations of eql3] we consider the inelastic 
scattering rate l/rin to be a constant independent on 
the electric field E and the electron energy e. 

Good agreement is found between the experiment and 
the numerical calculations for a broad range of temper- 
atures kT > r and magnetic fields. At small magnetic 
fields the conjecture of the independence of the inelastic 



time Tin on the electric field E is supported by direct 
evaluation of the variation (broadening) of the distribu- 
tion function, which is found to be small at the dc biases 
used in the experiment. The small variation provides a 
negligibly small correction to the inelastic collision in- 
tegral and to the inelastic scattering rate. Moreover at 
kT > r the energy space available for inelastic scatter- 
ing of an electron inside Landau sub-band contains, in 
fact, all levels of the sub-band. This may provide the 
weak dependence of the inelastic electron scattering on 
the energy e inside the Landau level. 

At a strong magnetic field, at which Landau levels are 
well separated, we have found a scaling behavior of the 
nonlinear resistance (see fig l9lf Op . In this regime the ex- 
periment and the theory demonstrate a remarkable cor- 
respondence even at a strong variation of the nonlinear 
resistance. This behavior is unexpected since the strong 
variation of the resistance implies a substantial deviation 
of the electron distribution function from the equilibrium 
and, therefore, an apparent inapplicability of the r ap- 
proximation with the constant Tin- Below we provide 
arguments, which shed a light on this interesting phe- 
nomenon. 

At a strong magnetic field, at which Landau levels are 
well separated, the spectral diffusion between Landau 
levels is absent due to the lack of the available electron 
states {f — 0). In this regime the total broadening of 
the distribution function is absent and, therefore, the to- 
tal number of Landau levels participating in the spectral 
diffusion is fixed. There is, however, a spectral diffu- 
sion inside Landau levels, generating local spectral flows. 
Since the spectral diffusion conserves the total number 
of particles and since there is no electron transport be- 
tween Landau levels, the total number of electrons inside 
any Landau level is preserved and equal to the thermal 
equilibrium value despite considerable deviations of the 
electron distribution function from the thermal equilib- 
rium inside the level. It is clear that in this condition 
the total number of empty states in each Landau level is 
also fixed and equal to the value at the thermal equilib- 
rium (at zero dc bias). Thus for the isolated Landau lev- 
els the averaged spectral distribution of electron states, 
which are available for the inelastic scattering of an elec- 
tron, is independent on the applied electric field. This 
may provide the significant stability of the inelastic relax- 
ation rate with respect to the dc bias. These arguments 
are valid, when the electron distribution inside a Landau 
level is not changing substantially with the electron en- 
ergy. This regime holds at relatively high temperature: 
kT > T. 

At low temperatures kT < T the only one Landau level 
is involved in electron transport and at the thermal equi- 
librium the electron distribution changes strongly inside 
the level. An application of a dc bias changes apprecia- 
bly the distribution of electrons. At kT < F the nu- 
merical calculations done in the r approximation deviate 
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substantially from the experiment (see figEJ;) , indicating 
a limited applicability of the approximation at the low 
temperatures. 

The numerical calculations are done in several steps. 
The goal of the first step is to find the density of elec- 
tron states i>{e) from a comparison with the experi- 
ment. The density of states ^(e) of the 2D electrons 
can be approximated by different theoretical expressions 
[H H [el [H, lil. We have found that the numeri- 
cal results for the temperature dependence of the inelas- 
tic scattering rate are robust with respect to particular 
choice of the expressions for the density of states (see 
below). Most of the numerical results, presented in the 
paper, are obtained using a Gaussian form of the DOS 
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iy{e) = i^Qy/uJ^'^ exp 



(e - nuJcY 

UJc/T^Tq 



(4) 



where the is the quantum scattering time. To find 
the DOS we compare normalized longitudinal resistance 
Rxxl Ro with the numerical evaluation of the normal- 
ized longitudinal conductivity Unil<JD obtained from eq[T] 
with thermal equilibrium distribution function fxi^)- 
The i?o is the resistance of the sample in zero magnetic 
field. In the leading approximation and at classically 
strong magnetic field {uJcTtr ^ 1) the two ratios equal to 
each other: R^x/Ro = o'ni/o'D- From the comparison we 
have obtained the quantum scattering time Tq and, there- 
fore, have approximated the density of electron states in 
eqlH Comparable values of quantum scattering time have 
been obtained using other methods, in particular, from 
analysis of magnitude of the quantum oscillations [58] . 

In the second step we use the DOS to numerically cal- 
culate the distribution function /(e) using eql3] in the 
limit t ^ Tin- In this limit the distribution function 
reaches a stationary state corresponding to the dc re- 
sponse. The distribution function is calculated at differ- 
ent values of the electric field E. 

In the third step the normalized nonlinear conductivity 
c^ni/cTD is calculated using eq[T]for different electric field. 
The results are compared with the normalized resistance 
Rxx/Rq- The inelastic scattering time is found from 
the best fit between dependencies of the normalized re- 
sistance Rxx/Rq and the calculated normalized conduc- 
tivity (Jni/aD on the dc bias. 

In accordance with eqO the spectral diffusion gener- 
ates an electron spectral flow from low energy re- 
gions (occupied levels) to high energies (empty levels). 
The spectral flow is proportional to the coefficient of the 
spectral diffusion and to the gradient of the distribu- 
tion function df /de: = D{e) ■ df /de. In a stationary 
state the spectral electron flow is constant. As a re- 
sult, the gradient of the distribution function df /de is 
strong in the regions of weak spectral diffusion (between 
Landau levels) and is small in the regions with strong 



FIG. 2: (color online) Normalized density of states u, distri- 
bution function / and non-equilibrium part of the distribu- 
tion function A/ = f — fr are shown as function of electron 
energy. The distribution function / is obtained by numeri- 
cal evaluation of eq. O using physical parameters typical for 
experiments presented below: Idc=377 {iJ,A); ri„=0.55 (ns); 
Tg^l.l (ps); B=0.924 (T) and T=10.7 (K) 



spectral diffusion (centers of the Landau levels). It is 
important to realize that a weak inelastic scattering can- 
not change significantly the robust dynamic flow in the 
energy space and, therefore, the behavior of the distri- 
bution function. This corresponds to our numerical cal- 
culations. Figl2] demonstrates the density of states, dis- 
tribution function and non-equilibrium part of the func- 
tion induced by dc current Idc- Indeed the gradient of 
the distribution function is considerably suppressed in- 
side Landau levels. This is due to both the fast spectral 
diffusion inside Landau levels and the slow diffusion be- 
tween them. Such non-equilibrium distribution function 
can not be described by a temperature 5J|. In accor- 
dance with eqU] the small gradient of the distribution 
function inside conducting Landau levels makes the net 
value of the nonlinear longitudinal conductivity (resistiv- 
ity) to be signiflcantly smaller than the linear, unbiased 
value. Below we present the detailed comparison between 
the experiments and the numerical calculations. 



RESULTS AND DISCUSSION 

Fig[3] demonstrates dependencies of the longitudinal 
resistance of two dimensional electrons on the magnetic 
fleld in sample N2. Two upper curves present dependen- 
cies obtained at different temperatures T=2.16K (dotted 
curve) and T=4.2K (solid curve) at zero dc bias. At small 
magnetic fields B <0.1T the magnetoresistance demon- 
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FIG. 3: (Color online), Dependencies of the longitudinal re- 
sistance Txx on magnetic field at different temperatures with 
no dc bias (black solid and dotted lines) and with applied dc 
bias Idc ~6 (mA) at T=2.04 K (grey solid line (red online)). 
Arrow indicates magnetic field B=0.1 T above which the elec- 
tron spectrum is modulated due to quantization of electron 
motion: Landau levels. 



strates the classical independence on the magnetic field 
[g^l- At B >0.1T the electron spectrum is quantized 
and at temperature T =0.3K the resistance demonstrates 
quantum oscillations (not shown). An arrow marks the 
magnetic field B =0.1T above which the electron spec- 
trum is modulated due to the quantization of the electron 
motion in magnetic fields. 

At magnetic fields B < 0.3T the two traces at 
T=2.16K and at T=4.2K are almost identical, indicat- 
ing a very weak temperature dependence of the resistance 
(drxx/dT > 0). At stronger magnetic fields the quantum 
oscillations (Shubnikov de Haas, SdH) are observed. The 
oscillations are result of Landau quantization of the elec- 
tron spectrum in the magnetic fields. At thermal equi- 
librium the amplitude A of the oscillations follows from 
eqUand eq[2]with the Fermi-Dirac distribution function: 
A ~ XT/sinh{XT), Xt = 2TT^kT/hoJc 



57, 581 . At smaU 



magnetic fields hujc <C kT the amplitude of the SdH os- 
cillations is small due to an effective averaging of the 
conductivity oscillations a{e) (see eqUJ over the temper- 
ature interval kT in eq[TJ Figl3] shows that the increase 
of the temperature reduces the magnitude of the oscil- 
lations symmetrically toward a background, which is an 
averaged value between maximums and minimums of the 
oscillations. 

A different behavior of the resistance is found in the re- 



sponse to the dc bias (5j| . In fig. 3 the lower curve presents 
a typical dependence of the differential resistance on mag- 
netic field at a finite dc bias. At B > O.IT, at which 



the Landau quantization appears, the resistance shows 
a considerable decrease with the dc bias [drxx/dl < 0). 
The decrease of the resistance cannot be explained by a 
temperature increase due to the dc heating. The temper- 
ature increase raises the resistance [drxx/dT > 0). More- 
over the quantum oscillations at the finite dc bias do not 
have the canonical shape, corresponding to the two up- 
per curves at zero dc bias. Instead a strong increase of 
higher harmonics of the oscillations is obvious. The en- 
hancement of the higher harmonic content is in apparent 
contradiction with the description of the dc biased elec- 
trons by an elevated temperature Tg: high temperature 
reduces exponentially the higher harmonic content of the 



54, 57 



oscillations 

Below we show that the strong decrease of the re- 
sistance with the dc bias is result of the non-uniform 
spectral diffusion of 2D electrons through Landau levels. 
We consider in detail two regimes. One regime corre- 
sponds to small magnetic fields, at which Landau levels 
are overlapped and the temperature is higher than the 
level separation: kT 3> fiujc- In this regime the quan- 
tum oscillations are absent and the resistance depends 
weakly on the temperature. At the small magnetic fields 
the spectral diffusion equation is solved both numerically 
and analytically '28'|. Another regime corresponds to high 
magnetic fields at which the Landau levels are separated: 
fiLOc > r. For sample N2 the first regime corresponds to 
B < 0.2T whereas the second regime is at S > 0.7T (see 
fig.©. 



Small magnetic fields 

At small magnetic fields the separation between Lan- 
dau levels hcuc is less than the effective width of the levels 
r = h/Tq. At low temperatures the width F is predom- 
inantly determined by the elastic impurity scattering of 
the 2D electrons. At small magnetic fields the density of 
states z/(e) is weakly oscillating with the energy e, mak- 
ing the spectral diffusion to also be a weakly modulated 
function of the energy. We consider a regime of high 
temperatures: kT 3> tiujc- In this regime the quantum 
oscillations are absent and the resistance increases weakly 
with the temperature T. 

Figdja) shows the dependence of normalized resistance 
R/Rq of the sample Nl on electric current at a small 
magnetic field B =0.343 (T) and temperature T =12.75 
(K). The parameter Rq is the resistance at zero mag- 
netic field. At small dc biases the normalized resistance 
decreases with the electric current. We consider the de- 
crease as a result of the non-uniform spectral diffusion 
of 2D electrons. At higher biases the resistance increases 
with the electric current due to other mechanisms of the 
nonlinearity 31 , S^l • In accordance with the theory (2^ 
the decrease of the resistivity obeys the following rela- 
tion: 
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T=12.75 (K); B=0.3434 (T); Sample 1 
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FIG. 4: (Color online), (a) Dependence of normalized lon- 
gitudinal resistance Rxx/{Ro ~ 37.75r2) on electric current . 
Symbols are experimental data points. Solid lines present an- 
alytical results (eqISj and numerical evaluation of the normal- 
ized resistance at 7 = 0.9931, = 1.138 (ps) and Tin =23.65 
(ps) for the gaussian form of the DOS. Thin dotted line is the 
numerical evaluation of the resistance, using the SCBA den- 
sity of states with 7 = 0.9931, = 1.132 (ps) and Tin =21.4 
(ps); (b) density of states, electron distribution function / 
and the non-equilibrium part of the function A/ = / — /t at 
dc bias Idc = 177-QfJ.A, (Gaussian DOS) ; (c) density of states, 
electron distribution function / and the non-equilibrium part 
of the function A/ = /-/t at dc bias he = 192. 5^, (SCBA 
DOS); T=12.75 (K), B=b.3434 (T), sample Nl. 



4Q 



dc 



i + Q, 



dc 



(5) 



where 7 = 1, (5 = exp{—'K/LOcTq) is the Dingle factor. The 
parameter Qdc takes into account the electric field E ( 
Hall electric field [H): 



2Tin 

Qdc = ( 

Ttr 



(6) 



To compare with the experiment we have used the Din- 
gle factor S{Tq) and the inelastic scattering time as 
fitting parameters. We also have varied parameter 7 to 
take into account possible memory eff'ects 27, [tO] and 
other deviations from the Drude magnetoconductivity 
7ll |. which are ignored at 7 = 1. A solid line presents the 
theoretical dependence (see eqlS)) of the normalized resis- 
tivity at 7 = 0.9931, Tq = 1.138 (ps) and Tin =23.65 (ps). 
Another solid line, which is indistinguishable from the an- 
alytical result, presents the numerical evaluation of the 
normalized resistivity, using eql3] with the same fitting 
parameters 7 = 0.9931, Tq = 1.138 (ps) and t^ =23.65 



(ps) and the Gaussian form of the DOS [6a|. A thin dot- 
ted line in figHJa) demonstrates the numerical evaluation 
of the resistance, using the SCBA density of states with 
7 = 0.9931, Tq = 1.132 (ps) and Tin =21.4 (ps). The 
density of states, electron distribution function / and 
the non-equilibrium part of the function A/ = / — /t 
are shown in figHb) (Gaussian DOS) and B^c) (SCBA 
DOS). FigHJa) demonstrates good agreement between 
the experiment and the theory at small dc biases. 

FigOJa) shows the dependence of the resistance of the 
sample N2 on the direct current at different tempera- 
tures as labeled. Solid lines present experimental depen- 
dencies. Dashed lines demonstrate results of numerical 
evaluation of the resistance, using eql3]with SCBA DOS 
at T=2.34 (K) and T=4.41 (K). The numerical calcu- 
lations demonstrate strong nonlinear suppression of the 
longitudinal resistance with the dc bias. The result is 
due to drastic modulation of the SCBA density of states 
and, therefore, spectral diffusion with the energy. 

The SCBA DOS, distribution function and the non- 
equilibrium part of the function are presented in the 
figEtb) at temperature T=4.41 (K). The DOS demon- 
strates sharp drops to almost zero values between Lan- 
dau levels. Such strong modulation of the DOS creates 
significant suppression of the energy exchange between 
different levels facilitating the electron " warming" inside 
the levels [H^l- The results, however, are apparently less 
compatible with the experiment than the one obtained 
with a smoother Gaussian DOS. 

In figlUJa) symbols present results of the numerical 
evaluation of the longitudinal resistivity, using eql3]with 
the Gaussian DOS and the quantum scattering times and 
inelastic times shown in the figEfa). The numerical sim- 
ulations demonstrate good agreement with the experi- 
ment in a considerably broader range of the dc biases. 
Gaussian DOS is shown in the figOJc), demonstrating 
moderate oscillations with energy. 

The experiment and the numerical calculations cor- 
respond well to each other at small electric currents 
Idc- At higher currents considerable deviations be- 
tween the experiment and the theory occur. The de- 
viations are expected. At higher currents there are 
additional mechanisms of the 2D electron nonlinearity 

[3, 0, [m, m 113, 0, m [ii, s, m, m, m, iia] , which 

are not taken into account in eqlH These nonlinearities 
are beyond the scope of the present paper. Moreover 
an additional contribution to the deviations may occur 
due to the conjecture of the constant inelastic relaxation 
rate 1/Ti„ in eqlS] At very small dc biases, at which the 
electron distribution is near the thermal equilibrium, the 
variation of the inelastic rate with the dc bias is also small 
since the phase space available for the inelastic scattering 
of an electron is nearly the same as at the equilibrium. At 
stronger dc biases the distribution function is broader and 
the inelastic scattering rate can be considerably stronger. 
To estimate the broadening of the distribution function 
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FIG. 5: (Color online), (a) Dependence of normalized lon- 
gitudinal resistance Rxx/Ro on electric current at diflterent 
temperatures as labeled. Solid lines are experimental curves. 
Symbols present result of numerical calculations of the resis- 
tance, using Gaussian DOS (eq|4} with 7=1 and and 
Tin presented in figl6ja); dotted lines demonstrate numerical 
evaluation of the R/Ro using SCBA DOS with 7=1 and Tq 
and Tin presented in fig|6ja). (b) Dependencies of normalized 
SCBA density of states i/(e) = //(e)/fo, electron distribution 
function / and non-equilibrium part of the function A/ on 
electron energy e counted with respect to Fermi energy /i. 
Distribution function a is solution of eq[3] using SCBA DOS 
with Tq=3.8 (ps), temperature T=4.41 (K) and electric cur- 
rent 7dc=50.6 (/xA). (c) Dependencies of normalized Gaussian 
density of states i>(e) — v{e)/vo, electron distribution function 
/ and non-equilibrium part of the function A/ on electron en- 
ergy e. The distribution function is a solution of. ([3] using the 
Gaussian DOS with rg=3.96 (ps), temperature T=4.41 (K) 
and electric current 1^^=56.4 (/iA); Ro{2.?,AK) = 44.6(r2), 
Roi'iAlK) = 46.36(n), Ro{6.17K) = 49.29(0), 7?o(8.41A') = 
52.47(0); B=0.2 (T); sample N2. 



at small magnetic fields, at which the spectrum is weakly 
modulated, we approximate the distribution function by 
an elevated temperature T^. At a stationary condition 
an increase of the Joule heat: dP = d{,P ■ p) is bal- 
anced by an increase of the heat dissipation: dE/Tr{Te) = 
c(Te)dT/Tr(Te), whcrc c(Te) = coTg is the electron heat 
capacity, is a time of the relaxation of the total elec- 



FIG. 6: Dependencies of the inelastic scattering time Tin and 
the quantum time Tq on temperature, (a) Filled squares show 
inelastic scattering time Tin, obtained numerically using eq|3] 
with Gaussian DOS; open circles present Tin obtained, using 
eqElwith SCBA DOS. Magnetic field B=0.2 (T). Sample N2. 
(b) Sample Nl. Gaussian DOS. Magnetic field is 0.5 (T). 



tron energy, J is current density and p is electron resis- 
tivity per square. In our case the time is controlled 
by the electron-phonon scattering, since the electron- 
electron scattering cannot stabilize the global broaden- 
ing of the distribution function. For the estimation of 
the broadening we use = Te-ph/T^ with Te-ph = 20 



(ns/K^) [60|, |77[. An integration of both sides of the 
balanced equation yields: — = 5Te-phJ^p/cQ. At 
the lattice temperature Tl=2.34 (K) the temperature in- 
crease ATe = Te-TL =0.14 (K) is found at Idc^9 (/lA). 
ATe=0.34 (K) is at Idc^lT (/iA), at which a deviation 
between the solution of cqOwith a constant Tin and the 
experiment is evident. Thus the estimation indicates that 
the deviation between the experiment and the theory at 
high dc biases can be also related to the variation of the 
inelastic scattering time Tin with the dc bias. Similar 
results are found for sample Nl. 

To obtain agreement between the experimental and 
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numerical dependencies in figl5^ we have used the con- 
stant inelastic scattering time Tin as a fitting parame- 
ter. The temperature dependence of the time Ti„, ob- 
tained from fitting at different temperatures, is shown 
in figlH for two samples. For sample N2 (figlUa) black 
squares) the inelastic time follows the dependence Tin — 
1.8(±0.3)/T2(±o i5) (ns). The time is obtained using 
Gaussian DOS shown in figElJc). Open circles in figlll^a) 
present the inelastic time Tin, obtained using the SCBA 
DOS shown in &gMh). The SCBA DOS results in con- 
sistently shorter inelastic times than the Gaussian DOS 
does, but with essentially the same temperature depen- 
dence. This holds for other magnetic fields and temper- 
atures. Taking into account the better overall agreement 
with the experiment obtained for numerical simulations 
with the Gaussian DOS, from now on we will only show 
numerical results for this density of states. 

Similar temperature dependence of the inelastic scat- 
tering time Tin is found for the sample Nl with a higher 
electron density and considerably shorter quantum scat- 
tering time Tq. The dependence is shown in ffglDJb). The 
dependence is obtained at magnetic field B=0.5 (T) and 
corresponds to the Gaussian DOS, which is similar to the 
one presented in figO^c). The quantum scattering times 
Tq in both samples are also shown for comparison and 
completeness in the figure. The time Tq is much shorter 
the inelastic scattering time Tin- The quantum scattering 
time has weak temperature dependence. 

In accordance with the theory the temperature de- 
pendence of the inelastic time ~ indicates the 
dominant contribution of the electron-electron scatter- 
ing into the inelastic relaxation of the distribution func- 
tion. We have compared the experimental results with 
theoretical calculations of the inelastic relaxation due to 
electron-electron interaction 2^, 7^ . For the param- 
eters corresponding to figE] the theoretical values of the 
inelastic time are found to be: t*^™^ = (ns) for 

sample N2 (figUa)) and r*^™'' = (ns) for sample 

Nl (figini^b)). The theoretical values are in good agree- 
ment with the experiment. A longer inelastic relaxation, 
found in the experiments, could be a result of an addi- 
tional screening by X-electrons in our samples [H^l- The 
screening is not taken into account in the comparison. 
Figl6] demonstrates a longer inelastic time for sample Nl 
with a higher electron density in agreement with the the- 
ory 



2i,[7i,l2: 



When considering the spectral diffusion of electrons in 
crossed electric and small magnetic fields at high tem- 
peratures, the results presented in this section demon- 
strate good quantitative agreement between the exper- 
iments and the theory. The numerical and analytical 
evaluation of the distribution function shows significant 
deviations of the electron distribution function from the 
Fermi-Dirac form leading to the nonlinear transport. At 
these conditions the rate of the inelastic relaxation of the 
non-equilibrium distribution function is found to be pro- 
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FIG. 7: (color online) (a) Relaxation of the non-equilibrium 
part of the distribution function A/ by an electron-electron 
scattering at small magnetic fields and/or high temperatures. 
Two electrons near maximum of A/ at energy eo scatter into 
nearest minimums at energies ei — eo — Ae and e2 = eo + Ae. 
The process conserves the total electron energy eo+eo = ei+e2 
and can be accomplished by the electron-electron interaction, 
(b) Inelastic relaxation at high magnetic fields and/or low 
temperatures. The relaxation fiows from overpopulated high 
energy levels (eo) toward under-populated low energy region 
(ei,e2). The relaxation flow does not conserve the total en- 
ergy of 2D electron system and cannot be accomplished by 
e — e scattering. The electron-phonon scattering provides the 
relaxation. 



portional to the square of the temperature: 1/Ti„ 



High magnetic fields 

At high magnetic fields the density of states and, there- 
fore, the spectral diffusion are strongly modulated with 
the energy. Between completely separated Landau levels 
(F <^ fiLOc) the spectral diffusion is expected to be very 
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weak. This may create a strong thermal isolation of the 
Landau levels and a stratification of the dynamic flow in 
the phase space in the response to the dc bias. In a lim- 
iting case of a single isolated level at low temperatures 
the global spectral flow is absent and the slope (gradient) 
of the distribution function df /de is determined solely by 
intra-level inelastic processes. For the intra-level inelastic 
transitions the electron-electron interaction may not be 
effective, because the interaction conserves the total en- 
ergy of electron system. FiglT] demonstrates a difference 
between the inelastic relaxation of distribution function 
through several Landau levels (figlUa)) and the relax- 
ation involving only one isolated Landau level (figlTlJb)). 

The first case (figlTja)) corresponds to a high tem- 
perature regime: kT 3> hujc- In the first case the 
electron-electron interaction can effectively reduce the 
non-equilibrium part of the distribution function A/ 
through the processes similar to the one shown in the 
figure. Two electrons near a maximum of the oscillating 
function A/ relax into the two nearest minimums. This 
process reduces the non-equilibrium part of the distri- 
bution function A/ smoothing out the oscillations. In 
this process the total electron energy is conserved and 
the relaxation can be accomplished by electron-electron 
scattering. 

The second case (figEl^b)) corresponds to low temper- 
atures (high magnetic field) fcT < F < huJc- Under 
these conditions the only Landau level (sub-band), lo- 
cated near the Fermi energy, is involved in the spectral 
diffusion. Lower energy levels are gapped and populated 
completely. They cannot participate in spectral trans- 
port due to the Pauli principle. The higher energy levels 
are empty, but, again, are inaccessible at low T due to 
the cyclotron gap. A typical non-equilibrium part of the 
distribution function corresponding to this case is shown 
in figElJb). The main flow of the relaxation to the ther- 
mal equilibrium is from overpopulated high energy levels 
into the under-populated low energy region of the Lan- 
dau level. The relaxation flow does not conserve the total 
energy of electron system, and, therefore, cannot be ac- 
complished by the electron-electron scattering. 

A possible candidate for inelastic electron relaxation is 
electron-phonon scattering. Electron-phonon scattering 
does not conserve the total electron energy and, there- 
fore, can be the mechanism responsible for the inelastic 
relaxation inside the isolated Landau level at low tem- 
peratures. Moreover, due to a stronger temperature de- 
pendence 7^ 77 1, the electron-phonon scattering could 
be the dominant mechanism of the relaxation at high 
temperature. Below we show the interplay between dif- 
ferent regimes of the inelastic electron relaxation, which 
are observed in our samples. 

FiglHUa) presents dependencies of the normalized re- 
sistance of the sample N2 at B = 0.784 (T) and at high 
temperatures as labeled. The magnetic field corresponds 
to a maximum of the SdH oscillations. At small currents 




FIG. 8: (color online) (a) Dependence of normalized re- 
sistance R/Ro on dc bias at high temperatures as labeled. 
Ro{6K) = 49.29(!:!), i?o(8.13is:) = 52.12(!:!). Insert demon- 
strates dependence of density of states, distribution function 
and non-equilibrium part of the function A/ on energy e; 
T=8.13 (K), he =58.5 (fiA), r„,=151 (ps),' =1.9 (ps). 
(b) Dependence of normalized resistance on dc bias at in- 
termediate temperatures from top to bottom at zero bias: 
T=1.48(iio = 43.68(n)), 1.97(iio = 44.33(0)), 2.44(i?o = 
44.99(0)), 2.93(7?o = 45.45(0)), 3.52(iio = 45.89(0)), 
4.08(7?o = 46.37(0)) (K). The electron system undergoes a 
transition to state with zero differential resistance at Idc > Ith 
and T <3 (K). Insert demonstrates dependence of density 
of states, distribution function and non-equilibrium part of 
the function Af on energy e; T=2.44 (K), he =18.2 (fiA), 
Ti„=3.77 (ns), Tq =2.75 (ps). (c) Dependence of normalized 
resistance on dc bias at low temperatures from top to bot- 
tom at zero bias: T=0.27(7io = 42(0)), 0.71(i?o = 42.64(0)), 
1.06(i?o = 42.99(0)) (K). Insert demonstrates dependence of 
density of states, distribution function and non-equilibrium 
part of the function Af on energy e; T=0.71 (K), he =6.67 
(/iA), ri„ = 17.7 (ns), Tq =3.65 (ps). Symbols are numerical 
calculations and solid lines are experiments. Magnetic field is 
0.784 (T). Sample N2. 
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the numerical simulation describes well the experiment. 
The insert to the figure shows the normalized density of 
states, distribution function / and non-equilibrium part 
of the function A/ at dc bias 58.5 (/iA). The regime cor- 
responds to the condition kT ^ F. 

FigHtb) presents dependencies of the normalized resis- 
tance at medium temperatures kT T. Again, at small 
currents the numerical simulation, obtained in the Tin ap- 
proximation of the right side of eqlH works well, provid- 
ing very good fit of the experiment data. At temperatures 
below 3 (K) a sudden deviation between the experimen- 
tal data and the simulation occurs above a threshold cur- 
rent of Ith= 6.6 (/iA). An arrow in the figure marks this 
current. It has been shown, that above the current Ith 
the electron system undergoes a transition into the zero 
differential resistance state 36|, |37| . In this state the dif- 
ferential resistance of the sample is nearly zero in a broad 
range of the current Idc > hh- Non- uniform, domain- like 
structures, propagating in real space, have been proposed 
to explain the origin of the electron state with zero dif- 
ferential resistance (27I [36j . Such states are beyond the 
regime described by the spatially uniform eqlS] 

It is interesting that the transition to the nonlinear 
state with zero differential resistance happens at a nor- 
malized value of the resistance Rtr = R/Ro ~ 1-5, which 
is almost independent on the temperature. Moreover 
at this point {Rtr j Ith) the nonlinear resistance demon- 
strates a transition from an insulating-like {dR/dT < 0) 
to a metallic-like {dR/dT > 0) behavior. These unex- 
pected features are currently not understood and will be 
subject of future studies. The insert to the figure shows 
the normalized density of states, distribution function / 
and non-equilibrium part of the function A/ obtained at 
dc bias 18.2 {fiA). 

Finally figlUJc) presents data at very low temperature 
kT < F. At this condition only one Landau level pro- 
vides the electron transport. At the low temperatures 
the theory, used in the Tin approximation, fits with the 
data only at very small currents. At the lowest temper- 
ature T=0.27K, numerical results deviate almost imme- 
diately from the experiment. The comparison indicates 
that the approximation of the inelastic collision integral 
in eqE] by a constant relaxation time Tin does not work 
in these conditions. At very low temperature the equilib- 
rium distribution changes very rapidly with the energy 
e inside the Landau level on a scale, which is much nar- 
rower than the level width F: kT <C F. Since the inelastic 
processes are extremely weak at the low T, the spectral 
diffusion broadens easily the electron distribution to a 
scale comparable with the width of the level F even at 
small dc biases. This process increases significantly the 
phase space available for the inelastic electron scattering, 
enhancing the scattering rate 1 /Tin appreciably. Thus at 
kT < F the inelastic scattering depends strongly on the 
dc bias and the spectral diffusion equation (eq[3]) with a 
constant ri„ does not describe the nonlinear resistance 



appropriately. More work is required to evaluate quan- 
titatively the shape of the distribution function in this 
regime. However we suggest that even in the regime 
kT < F the distribution function will be qualitatively 
similar to the one shown in the insert to figEI^c), which 
is obtained in the r approximation. At a high dc bias 
the function can not be described by an elevated elec- 
tron temperature as it is shown in the figure (see also 

Additional analysis of the curves at the high magnetic 
fields reveals an interesting scaling behavior of the non- 
linear resistance. Applying two linear transformations 
{y — Ky ■ y and x = • x) along y and x-axes one 
can collapse all dependencies at different temperatures 
presented in figl51[a,b) on a single curve. FiglH^a) shows 
the result. The y-transformation normalizes the resis- 
tance at zero bias to unity: R{I) = R{I)/R{I = 0). The 
linear x-transformation, applied along the x-axes, pro- 
vides the final result. Solid curves are experimental de- 
pendencies measured in temperature interval (1.48-8.13) 
(K). Open circles show a result of numerical calculations 
of the nonlinear resistance obtained using eql3] with the 
equilibrium electron distribution at T =4.08 (K) and 
Tq = 2.75 (ps). The same scaling is found for sample 
Nl in a broader range of temperatures. The result is 
shown in fig fTUT a). All dependencies are plotted versus 
a parameter A^/"^ — {a^^E'^Tin/ voY^"^ ^ Idc- At a fixed 
density of states v{e) the variable A ~ E^Tin is the main 
parameter, which determines the deviation of the elec- 
tron distribution / from the thermal equilibrium fx in 
eqU 

FiglQlJa) demonstrates a good scaling and a remark- 
able correspondence with numerical results obtained at 
A^/"^ < 0.15, using eql3] with a fixed Tin- The cor- 
respondence between the experiment and the theory is 
even more impressive for a curve at the lowest tempera- 
ture (T=2.34 (K)) presented in figflWa). Almost perfect 
agreement between the experiment at T=2.34 (K) and 
the theory is found at substantially stronger dc biases 
{A ~ 1). The scaling of the nonlinear resistance and the 
excellent agreement with the theory indicates strongly 
the presence of the spectral diffusion with a constant rate 
of the inelastic relaxation l/ri„. 

We suggest that the scaling is a result of a specific non- 
linear regime, which occurs for separated Landau levels. 
As we have already mentioned in the section " Theory and 
Numerical Simulations" , the spectral diffusion between 
well-separated Landau levels is absent. In this regime 
there is no global broadening of the distribution function. 
Moreover inside each of the Landau levels the local spec- 
tral flow preserves the number of electrons and, there- 
fore, the number of the empty states. Thus the strati- 
fied spectral diffusion keeps the spectral distribution of 
the available phase space (averaged over each Landau 
level), to be fixed and the same as the one at the ther- 
mal equilibrium {E — 0). The invariance of the phase 
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FIG. 9: (a) Scaling of normalized resistance with parameter 
A°'^ ~ /dc- AH curves presented in figl8ja,b) at diflterent tem- 
peratures (1.48-8.13) (K) follow the same dependence on the 
parameter yl"'^ < 0.15 (solid curves). Open circles present 
results of numerical calculations of the normalized resistance, 
using eqOJwith = 2.75 (ps), T = 4.08 (K), B=0.784 (T) 
and parameter A^^^ = (ajj'^E^Tin/i/oy^^', insert shows inde- 
pendence of variations of the normalized resistance with A 
on temperature T. The results are obtained using eq|3] at 
T=3(K) -open circles, T=4.08(K) - solid curve, and T=6(K)- 
fiUed circles, (b) Dependences of inelastic scattering time Tin, 
obtained from comparison between experiment and numerical 
evaluation of nonlinear resistance, using eq[3] (filled squares) 
and from scaling (open circles) on temperature. Open squares 
present temperature dependence of quantum scattering time 
Tq. Magnetic field 5=0.784 (T). Sample N2. 



space available for inelastic processes could provide the 
independence of the inelastic scattering time rj„ on the 
dc bias fixing the time at the thermal equilibrium value: 
Tin{E) = Tin{E = 0). The constant inelastic scattering 
rate makes the evolution of the electron distribution and 
the nonlinear resistance to be universal in a broad range 
of the dc biases. 

The scaling reveals another interesting property of the 



FIG. 10: (a) Scaling of normalized resistance (solid curves) 
with parameter A'^'"' ~ Idc at different temperature from bot- 
tom to top: 2.34, 4.2, 5.4, 7.6, 10.7, 14.8, 20.1, 24.6 (K). 
Open circles present results of numerical calculations of the 
normalized resistance, using eq|3]with Tq — 1.1 (ps), T — 2.34 
(K), 5=0.924 (T) and parameter A'^^'^ = (cr£5Vi„/i/o)^''^; 
(b)Dependences of inelastic scattering time Tin, obtained from 
comparison between experiment and numerical evaluation of 
nonlinear resistance, using eq|3] (filled squares) and from scal- 
ing (open circles) on temperature. Open squares present tem- 
perature dependence of quantum scattering time Tq. Magnetic 
field 5=0.924 (T). Sample Nl. 



nonlinear regime. FiglH^a) shows that variations of the 
normalized resistance with parameter A^^'^ <0.15 is the 
same at different temperatures and, therefore, does not 
depend on the initial, equilibrium distribution fx of 2D 
electrons in eqlJ] The equilibrium distribution fx is sub- 
stantially different in the temperature interval, in which 
the scaled dependencies have been measured: (1.4 - 8.13) 
(K). We suggest that the independence of the nonlin- 
ear resistance on the /t is also a result of the absence 
of the dc bias induced spectral flows between Landau 
levels. Without the inter-level spectral flow the levels 
are, in essence, independent from each other and, there- 
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fore, absorb the energy from electric field independently. 
The absorption inside each Landau level is determined 
by the same spectral dynamics, assuming that the den- 
sity of states is the same for each level. An estimation 
of the nonlinear conductivity in a model of separated 
(independent) levels supports the suggestion [7J|. The 
numerical evaluation of the nonlinear behavior of the re- 
sistance, which has been done for different temperatures, 
using eql3l demonstrates also the independence of the 
normalized nonlinear resistance on the temperature in 
this regime. In particular, the numerical values of the 
normalized resistance obtained for T=3K, T=4.08K and 
T=6K at a fixed density of states (r^ = 2.75 (ps)) differ 
by less that 3% at any A <0.4. This is shown in the 
insert to figElJa). 

The scaling of the nonlinear resistance provides an easy 
practical access to the variation of the inelastic relax- 
ation time with the temperature since it does not re- 
quire the solution of the eql31 The scaling coefficient 
Kx ~ E ■ (Ti„(T))^/^ takes into account the temperature 
variations. A comparison of the inelastic time Tin ob- 
tained from the scaling (open circles) and from the direct 
comparison with the numerical calculation of the nonlin- 
ear resistance using eq|3] (solid squares) are presented in 
figlUb) (sample N2) and figdUIb) (sampleNl). There is 
a good overall agreement between two approaches. A dif- 
ference appears since the numerical calculation takes into 
account a variation of spectral dynamics with the tem- 
perature due to changes in density of states (see the time 
Tq presented in the figures) and a temperature variation 
of the transport scattering rate. 

Deviations from the scaling depend on the tempera- 
ture. Presented in figlUa) and fig fTUf a) at higher tem- 
peratures experimental curves deviate up from the scal- 
ing behavior at a smaller A. Taking into account the 
strong reduction of the inelastic scattering time with 
the temperature, one can find that the deviations from 
the scaling occur at progressively higher dc biases: E ~ 



{A/, 



U/2 



This indicates that corrections to the scaling 



due to other nonlinear mechanisms, arising at high biases 
0, mi nil , decreases with the temperature increase. The 
later agrees with the temperature dumping of a magni- 
tude of the dc bias induced magneto-oscillations of the 
nonlinear resistance due to inter- level scattering 
At high dc biases A^^"^ > 0.15 sample N2 demonstrate 
an additional abrupt deviation down from the scaling at 
temperatures below 3K (see figUJa)). As we have men- 
tioned at this condition a transition to the zero differen- 
tial resistance state appears (36l [37j . which may break 
down the description of the 2D electron system by the 
spatially uniform spectral equation (eq[3]) |36j . 

Below we discuss the temperature dependence of the 
inelastic scattering time. Figinb presents the temper- 
ature dependence of the time at magnetic field 
5=0.784 (T) for the sample N2. Two temperature 
regimes are clearly observable. At temperatures T > 2K 



the inelastic relaxation time is inversely proportional 
to T^: Tin = 66(±10)/T3(±o i5) (ns). At temperatures 
below 2K the inelastic time depends weaker on the tem- 
perature: r,„ = 11.6(±2)/Ti-26±o i5 (ns). 

The observed dependence of the inelastic time 
Tin = 66/r^(ns) correlates with the one obtained in Si- 
MOSFETs : n„ = (10 - 60) /T^ (ns) at temperatures 
1.5 < T < 4.2K [Hi] and with the dependence found 
in a GaAs/AlGaAs heterojunction: Tin — 2G/T'^(ns) 
at temperatures I < T < 3K [g^I- In both papers 
the temperature dependence has been attributed to the 
electron-phonon scattering. We suggest that the tem- 
perature dependence observed at T >2K is also due to 
an electron-phonon scattering in Bloch-Gruneisen (BG) 
regime at which the wave vector of a typical thermal 
phonon qx = kT/hs is smaller than the size of the Fermi 
circle 2k p: qt < 2k p. Here s is sound velocity and kp 
is Fermi wave vector (63| . In our high density samples 
the BG regime exists at temperatures below Tbq « 20K, 



where kTsc — 2fci? • fis 75[. A theoretical evaluation 



of the inelastic electron-phonon scattering time in GaAs 
quantum wells due to screened piezoelectric (PZ) cou- 
pling yields: tpz ~ 1&/T^ (ns) at temperatures of few 
K at zero magnetic field [tR [tH . Deformation potential 
(DP) yields a comparable contribution to the electron- 
phonon scattering rate at T >4K. At a weak screen- 
ing the electron-phonon scattering time is found to be 
Top ^ 18/T^(ns) [ill at zero magnetic field. 

The temperature dependence is found also for 

the sample Nl at high temperatures. FigllOb presents 
the temperature dependence. At T >10K the inelas- 
tic scattering time is proportional to IjT^: Tin = 
70(±10)/T3=^°-2 (ns). The dependence is the same as 
the one observed in the sample N2. At lower tempera- 
tures T < 10 (K) the inelastic relaxation time deviates 
consistently from the dependence. The temperature 
dependence Tin — 9(±2)/T^(='=°-^) provides a reasonable 
approximation, indicating a possible contribution of the 
electron-electron interaction to the inelastic relaxation 
rate. The same (T"^) temperature dependence is ob- 
served at small magnetic fields for both samples but at 
considerably stronger relaxation rate. Thus the tempera- 
ture dependence below 10(K) appears as an intermediate 
regime at which the electron-electron scattering is signif- 
icant but is suppressed considerably by the quantization 
of the electron spectrum. At the beginning of the section 
we have discussed the possible reason for the reduction 
of the contribution of the e — e scattering to the inelastic 
relaxation in strong magnetic fields. 

Our experiment demonstrates a correlation between 
modulation of the density of states, the inelastic time 
Tin and the temperature dependence of the time. At 
low magnetic field B=0.2 (T) the density of states of 
the sample N2 is weakly modulated at about ±40% (see 
figO. The time of inelastic relaxation equals to 1.8/T^ 
(ns) below 8K. At the magnetic field B =0.784 (T) the 
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FIG. 11: Dependences of inelastic scattering time nn and 
quantum scattering time r, on magnetic field at two different 
temperatures as labeled. Two shaded areas indicate two dif- 
ferent temperature regimes of the inelastic electron relaxation 
observed in the sample. Sample N2. 

modulation of the density of states of the sample N2 is 
significantly stronger approaching 95 % of the averaged 
value (see figlSl[a,b)). The inelastic time equals to 66 /T^ 
at 2 < T < 8 (K). In magnetic field B =0.924 (T) the 
modulation of the density of states of the sample Nl is 
about 60% and the inelastic time is between the two pre- 
vious values: l.S/T^ <9/T^ < 66 /T^ at T < 7 (K). 

In accordance with the correlation one should expect a 
gradual reduction of the contribution of electron-electron 
scattering to the inelastic relaxation and an increase of 
the relaxation time rj„ with an increase of the modula- 
tion of the density of states. An increase of the magnetic 
field B enhances the DOS modulation. FigHT] presents 
the dependence of the inelastic time Tin on the magnetic 
field for sample N2 at two different temperatures as la- 
beled. Magnetic field increases the relaxation time Tm- 
The temperature dependence of the inelastic relaxation 
rate changes from at low magnetic field to at high 
magnetic fields. In the figure, two rectangular shaded ar- 
eas indicate the two different temperature regimes of the 
inelastic relaxation. These regimes are presented in more 
details in figElJa) and figlljb). Similar enhancement of 
the relaxation time Ti„ with the increase of the magnetic 
field is found for sample Nl (not shown). 

CONCLUSION 

We have studied the nonlinear response of 2D electrons 
placed in crossed electric and quantized magnetic fields 
at low temperatures. The resistance of 2D electrons de- 
creases strongly with an increase of the electric field. The 
decrease of the resistance is in good quantitative agree- 



ment with theory considering the nonlinear response as 
a result of non-uniform spectral diffusion of 2D electrons 
limited by inelastic electron scattering. Comparison be- 
tween the experiments and the theory has revealed dif- 
ferent regimes of the electron inelastic relaxation. 

At low magnetic fields, at which the Landau levels are 
well overlapped and the spectral diffusion is weakly mod- 
ulated with the electron energy, the inelastic scattering 
rate is found to be proportional to the square of the tem- 
perature in temperature interval (2-10 (K)). The de- 
pendence indicates the electron-electron scattering as the 
dominant mechanism of the inelastic relaxation. At high 
magnetic fields, at which the Landau levels are well sep- 
arated, the spectral diffusion is strongly modulated and 
the rate of the inelastic relaxation is proportional to T^. 
This suggests the electroii-phonon scattering to be the 
dominant inelastic mechanism. At fixed temperature the 
inelastic time Tin increases with the magnetic field. At 
very small temperatures kT < F and well separated Lan- 
dau levels an additional regime of the inelastic electron 
relaxation is identified: 1/Ti„ ~ T^'^^. 

At the high magnetic fields the nonlinear resistance 
demonstrates scaling behavior in a broad range of tem- 
peratures exceeding the width of Landau levels. The scal- 
ing indicates specific regime of the dc heating in electron 
systems with discrete electron spectrum. A temperature 
cannot describe the heating. The spectral diffusion lim- 
ited by the inelastic relaxation with constant rate de- 
scribes remarkably well the scaling in broad range of the 
dc biases. 
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